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Abstract: In this short survey we recollect some of the recent results on the high 
energy behavior (i.e., for diverging sequences of eigenvalues) of nonlinear function¬ 
als of Gaussian eigenfunctions on the d-dimensional sphere d > 2. We present a 
quantitative Central Limit Theorem for a class of functionals whose Hermite rank 
is two, which includes in particular the empirical measure of excursion sets in the 
non-nodal case. Concerning the nodal case, we recall a CLT result for the defect 
on . The key tools are both, the asymptotic analysis of moments of all order for 
Gegenbauer polynomials, and so-called Fourth-Moment theorems. 
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1 Introduction 

Let us consider a compact Riemannian manifold (A4,g) and denote by its 
Laplace-Beltrami operator. There exists a sequence of eigenfunctions {/j}jeN and 
a corresponding non-decreasing sequence of eigenvalues 

fj + Ejfj = 0 ) 

such that {fj}j£N is a complete orthonormal basis of the space of square 

integrable measurable functions on M. One is interested in the high energy be¬ 
havior i.e., as j —>■ -|-oo, of eigenfunctions fj, related to the geometry of both level 
sets fj~^{z) for z G R, and connected components of their complement 
One can investigate e.g. the Riemannian volume of these domains: a quantity that 
can be formally written as a nonlinear functional of fj. 

The nodal case corresponding to z = 0 has received great attention (for moti¬ 
vating details see [11]). 

At least for “generic” chaotic surfaces A4, Berry’s Random Wave Model allows 
to compare the eigenfunction fj to a “typical” instance of an isotropic, monochro¬ 
matic random wave with wavenumber (see [11]). In view of this, much effort 
has been first devoted to 2-dimensional manifolds such as the torus (see e.g. 
[4]) and the sphere (see e.g. [3], [2], [8], [12]). Spherical random fields have 
attracted a growing interest, as they model several data sets in Astrophysics and 
Gosmology, e.g. on Cosmic Microwave Background ([5]). 

More recently random eigenfunctions on higher dimensional manifolds have 
been investigated: e.g. on the hyperspheres ([6]). 
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1.1 Random eigenfunctions on 

Let us fix some probability space (11, P), denote by E the corresponding expec¬ 
tation and by C the unit d-dimensional sphere (d > 2); stands for the 

Lebeasgue measure of the hyperspherical surface. By real random field on we 
mean a real-valued measurable map defined on (17 x (8) 7?(S‘^)), where 

denotes the Borel cr-field on Recall that the eigenvalues of the Laplace-Beltrami 

operator A§d on are integers of the form —£{£ -|- d — 1) =: —Ei, t' G N. 

The f'-th random eigenfunction Tg on is the (unique) centered, isotropic real 
Gaussian field on S'* with covariance function 


Ke{x,y) := Gi.d{cosT{x,y)) a;,yGS‘*, 


where Gt^d stands for the f-th Gegenbauer polynomial normalized in such a way 
that G^;d(l) = 1 and r is the usual geodesic distance. More precisely, setting 


at,d 


:= ( +1 ^), we have G^-d = 


’-e-d 


(f-Lf-i) 


where P} 


»./ 3 ) 


denote standard 


Jacobi polynomials. By isotropy (see e.g. [5]) we mean that for every g G S'O(d-l-l), 
the random fields = {Ti{x))x^‘^d and T® := (Ti(gx))^^§d have the same law in 
the sense of finite-dimensional distributions. Here SO{d + 1) denotes the group 
of real (d -I- 1) x (d -I- l)-matrices A such that AA' = I the identity matrix and 
detH = 1. 


Random eigenfunctions naturally arise as they are the Fourier components of 
those isotropic random fields on S'* whose sample paths belong to L^(S'*). 

Let us consider now functionals of Ti of the form 


St{M) := f M{Tt{x))dx , 


( 1 ) 


where M : R ^ R is some measurable function such that E[M(Z)^] < -|-oo, 
Z ~ A/’(0,1) a standard Gaussian r.v. In particular, if M(-) = !(• > z) is the 
indicator function of the interval (z,-|-oo) for z G R, then (1) coincides with the 
empirical measure Si{z) of the z-excursion set Ai{z) := {cc G S'* : Ti{x) > z}. 


1.2 Aim of the survey 

We first present a quantitative GLT as f —>■ -l-oo for nonlinear functionals Si{M) in 
(1) on S'*, d > 2, under the assumption that ¥\M{Z)H 2 {Z)] ^ 0, where H 2 {t) := 
7^ — 1 is the second Hermite polynomial. 

For instance the above condition is fullfilled by the empirical measure Si{z) of 
z-excursion sets for z yf 0. For the nodal case which corresponds to the defect 

De := [ liTeix) > 0) dx - [ l{Te{x) < 0) dx , ( 2 ) 

we present a GLT for d = 2. Quantitative CLTs for on S'*, d > 2, will be treated 
in a forthcoming paper. 

We refer to [7], [8] and [6] for the spherical case d = 2 and to [6] for all higher 
dimensions. The mentioned results rely on both, the asymptotic analysis of mo¬ 
ments of all order for Gegenbauer polynomials, and Fourth-Moment theorems (see 

[9], [1]). 
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2 High energy behavior via chaos expansions 

For a function M : K. —s- R as in (1), the r.v. Si{M) admits the chaotic expansion 

SiiM) = V f H,iTi{x)) dx (3) 

(see [9]) in L^(P) (the space of finite-variance r.v.’s), where Hq is the g-th Hermite 
polynomial (see e.g. [10]) and Jq{M) := ¥,[M{Z)Hq{Z)\, Z ^ A/’(0,1). We have 
E[Se{M)] = Jo{M)fXd] w.l.o.g. Jo(M) = 0. 

The main idea is first to investigate the asymptotic behavior of each chaotic 
projection, i.e. of each (centered) r.v. of the form 

hi-q,d-= [ Hq{Ti{x))dx (4) 

and then deduce the asymptotic behavior of the whole series (3). Note that he-i^d = 
0, as Tf has zero mean on By the symmetry property of Gegenbauer polynomials 
([10]), from now on we can restrict ourselves to even multiples 1', for which some 
straightforward computations yield 

^7r/2 

Var[hi.q^d] = 2ql^ld^^d-l G^;d(cosd)‘^(sind)‘^“^ dd . (5) 

do 


2.1 Asymptotics for moments of Gegenbauer polynomials 


The proof of the following is in [7], [8] for d = 2 and in [6] for d > 3. 
Proposition 1. As £ ^ oo, for d = 2 and q = 3 or q > 5 and for d,q >3, 


Gf;d(cosd)«(sind)'*-' dd = (1 + o(l)) . 


The constants Cq-^d are given by the formula 


( 6 ) 


^q;d • — 






+d-l 


dif , 


(7) 


where Jd_^ is the Bessel function ([10]) of order f — 1. The r.h.s. integral in (7) is 
absolutely convergent for any pair (d, q) (2, 3), (3,3) and conditionally convergent 
for d = 2,q = 3 and d = q = 3. Moreover for C 4;2 := 


Gf; 2 (cosd)"‘sinddd = C 4 ; 2 ^^^(l-b o(l)) . (8) 

From [10], as ^ -boo, 

J G^;d(cosd)^(sind)‘^“^dd = 4/XdAid-i'^|i^(l-b o(l)) , C 2 ;d := • (9) 

Clearly for any d, q > 2, Cq-d > 0 and Cq-d > 0 for all even q. Moreover we can give 
explicit expressions for C 3 ; 2 ,C 4;2 and C 2 -d for any d > 2. We conjecture that the 
above strict inequality holds for every pair (d, g), and leave this issue as an open 
question for future research. 
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2.2 Fourth-Moment Theorems for chaotic projections 

Let us recall the usual Kolmogorov dx, total variation drv and Wasserstein dw 
distances between r.v.’s X,Y: for G {K, TV,W} 


dv{X,Y)-.= sup |E[/r(X)]-E[/i(y)]| , 

hGE-D 


where Hk = {!(• < z),z € M}, Htv = {1a(-)M G B(M)} and Hw is the set of 
Lipschitz functions with Lipschitz constant one. 

The r.v. he-q^d in (4) belongs to the so-called gth Wiener chaos. The Fourth- 
Moment Theorem ([9]) states that if Z ~ A/’(0,1), for V G {K, TV,W} we have 






^/ya.T[hi.q^d 


■ ,z <Cc(g)i 


I cuin4{he.q^d) 

Var[h^;,,d]2 


( 10 ) 


where C-piq) > 0 is some explicit constant and cwo[i 4 {ht,q,d) is the fourth cumulant 
of the r.v. hi-q^d- An application of (10) together with upper bounds for cumulants 
leads to the following result (see [6]). 

Theorem 1. For all d,q>2 and V G {K, TV, W} we have, as £ ^ -foo, 

dr, ( ,z] = O (log£)-"(«’'')) , (11) 

\^yXar[hi.q^d] J ^ ^ 

where d{q; d) G Q, r]{q; d) G {—1, 0,1} and r]{q; d) = 0 but for d = 2 and q = 4,5, 6. 

The exponents S{q;d) and r]{q;d) can be given explicitly (see [6]), turning out 
in particular that if {d, q) ^ (3, 3), (3,4), (4, 3), (5, 3) and Cq-d > 0, 

he.,q,d - ^ Z , as £ ^ -foo , (12) 

^yYar[hi.q^d] 

where from now on, denotes convergence in distribution and Z ~ A/"(0,1). 

Remark 1. For d = 2, the CLT (12) was already proved in [8]; nevertheless Theo¬ 
rem 1 improves the existing bounds on the rate of convergence to the asymptotic 
Gaussian distribution. 


2.3 Quantitative CLTs for Hermite rank 2 fnnctionals 

Proposition 1 states that whenever M is such that J 2 (M) ^ 0 in (3), i.e. the 
functional Si{M) in (1) has Hermite rank two, then 


lim 

£—¥-\-00 


Var[Se{M)] 


Var 


J2(M) 

2 


hl-2,d 


= 1 . 


(13) 


Hence, loosely speaking, S({M) and its 2nd chaotic projection have 

the same high energy behaviour. The main result presented in this survey is the 
following, whose proof is given in [6]. 
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Theorem 2. Let M : R R in (1) be s.t. E [M{Z)H 2 {Z)] =: J 2 {M) ^ 0, then 

dwj - ,z\=o(£-i] , asi^oo, (14) 

\^VaT[S,{M)y J \ J ’ 

where Z ~ A/’(0,1). In particular, as i ^ +oo, 

4 Z . (15) 

VVar[^,(M)] 


3 Geometry of high energy excursion sets 


Consider the empirical measure Se{z) of the z-excursion set Ai{z) for z G R, as 
in §1.1. It is easy to check that in (3) E[5'£(z)] = ^^(1 — <i>(z)) and for g > 1, 
Jg(l(- > z)) = Hq-i{z)4>{z), where $ and </> denote respectively the cdf and the pdf 
of the standard Gaussian law. Since J 2 (l(' > z) = z(f>{z), Theorem 2 immediately 
entails that, as £ ^ oo, if z 0 


dw 


I Se{z)- Hdjl - ^{z)) 
\ \/Var[S'^(z)] 



The nodal case z = 0 requires different arguments: in the chaos expansion for the 
defect (2) Di only odd chaoses occur but each of them “contributes” by Proposition 
1. Asymptotics for the defect variance on have been given in [7]: 

C 

Yai[Di] = — (1 + o(l)) , as £ -)■ +oo , 
for C > ■^=. Moreover in [8] a CLT has been proved: as £ —k +oo, 


4 Z 

i/Varp^ 

where Z ~ ^^(0,1). In a forthcoming paper, we will provide quantitative CLTs for 
the defect on S'^, d > 2. 

Remark 2. The volume of excursion sets is just one instance of Lipschtz-Killing 
curvatures. In the 2-dimensional case, these are completed by the Euler-Poincare 
characteristic ([3]) and the length of level curves ([4],[12] for the nodal variances). 
In forthcoming papers jointly with D. Marinucci, G. Peccati and I. Wigman, we 
will investigate the asympotic distribution of the latter on both the sphere and 
the 2-torus T^. For future research, we would like to characterize the high energy 
behavior of all Lipschitz-Killing curvatures on every “nice” compact manifold. 
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